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We examine the consequences of the existence of correlated mixtures of adiabatic and isocurvature 
perturbations on the CMB and large scale structure. In particular, we consider the four types of 
"elementary" totally correlated hybrid initial conditions, where only one of the four matter species 
(photons, baryons, neutrinos, CDM) deviates from adiabaticity. We then study the height and 
position of the acoustic peaks with respect to the large angular scale plateau as a function of the 
isocurvature to adiabatic ratio. 



I. INTRODUCTION 



In the near future, a lot of data about the anisotropies of the Cosmic Microwave Background (CMB) will be available 
to cosmologists, notably thanks to balloon experiments and the planned satellites MAP [jlj and Planck j2|. What will 
be remarkable is the expected high resolution and sensitivity of these experiments, which may turn cosmology into a 
high precision activity. 

One of the hopes of cosmologists is to be able to determine from these data the cosmological parameters describing 
■ the geometry and matter contents of our Universe. In this respect, it is important to stress that the fluctuations that 
are and that will be measured, result, according to our current understanding, from a combination of the primordial 
0^ ■ perturbations and of the cosmological parameters. In the preparation of the future data analysis, one should be careful 
to avoid oversimplification a priori of the primordial perturbations and not to stick to the simplest one-scalar field 
inflation model. After all, the early universe is the period in the history of the Universe where the physics is the least 
' known. 

, A more general description of the primordial perturbations may therefore be needed to be able to interpret the 
future data. In this perspective, the aim of this work is to examine the consequences of the existence of isocurvature 
perturbations in addition to the usual adiabatic perturbations. Such studies have already been performed in the case 
of independent mixtures of adiabatic and isocurvature perturbations ||- S. This is why we will focus our attention 
on correlated mixtures of adiabatic and isocurvature perturbations. The possibility of such primordial perturbations 
+3 , is motivated by the recent work of one of us, which showed that the simplest model of multiple inflation, a model 
with two massive non-interacting scalar fields, can produce such correlated mixtures [fTo) . 

Isocurvature perturbations are perturbations in the relative density ratio between various species in the early 
universe, in contrast with the more standard adiabatic (or isentropic) perturbations which are perturbations in the 
total energy density with fixed particle number ratios. Primordial isocurvature perturbations are often ignored in 
inflationary models. The main reason for this is that they are less universal than adiabatic perturbations because, 
on one hand, they can be produced only in multiple inflationary models and, on the other hand, they do not 
necessarily survive until the present epoch. 

However, isocurvature perturbations have been shown to be of potential importance in some specific models : 
axions @, §], 0, |l§], Affleck-Dine baryogenesis mechanism (n), multiple field inflation @, [|fl), Q, 



A priori, since the (not too old) Universe is filled with four species, baryons, photons, neutrinos and dark matter 
(which will be assumed to be cold here), several types of isocurvature perturbations can be envisaged. For example, 
in the past, a model with isocurvature baryon perturbations was proposed fl(i|| , although it does not seem compat- 
ible with the data today ]l7j . Most of the recent models, however, contain cold dark matter (CDM) isocurvature 
perturbations. A more general approach, including neutrino isocurvature perturbations (and also isocurvature veloc- 
ity perturbations), was considered recently plf . In the present work, we will focus our attention, for simplicity, on 
primordial perturbations where only one species deviates from adiabaticity, which thus leaves room to four types of 
hybrid (i.e. adiabatic plus isocurvature) initial perturbations. These four "elementary" modes will be systematically 
studied, without trying to make any connection with specific early universe models. 

As far as observational constraints are concerned, it has already been established that a pure isocurvature scale- 
invariant spectrum must be rejected because it predicts on large scales too large temperature anisotropies with 
respect to density fluctuations but other possibilities have been envisaged, like tilted isocurvature perturbations. 
The main trend, however, has been to study models with a mixing of isocurvature and adiabatic perturbations. 
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Confrontation of these models with observational data, like CMB anisotropics and large scale structures, seems to 
allow only for a small fraction of isocurvature perturbations. Future CMB measurements will also enable us to put 
much tighter constraints on this kind of models. 

It must be emphasized, however, that all these studies assumed independent mixing of isocurvature and adiabatic 
perturbations. While this assumption can be indeed justified in some specific early universe models, it is certainly not 
an absolute rule, as it has been shown in |ic[ |. It is thus the purpose of this paper to investigate the consequences on 
observational quantities, namely the large scale structure and the CMB anisotropics, of correlated mixtures of isocur- 
vature and adiabatic perturbations. As it will be shown, correlation gives more richness to hybrid perturbations. For 
example, while the first acoustic peak (relatively to the plateau) is always lower for independent hybrid perturbations 
than for pure adiabatic perturbations, it can be either lower or higher with correlated hybrid perturbations. 

The plan of the paper will be the following. In the next section (§||), we recall the basic definitions of isocurvature 



and adiabatic perturbations and introduce hybrid perturbations. Then, §111 will discuss the notion of correlation 



between isocurvature and adiabatic perturbations. In { IV, we will begin the systematic analysis of the correlated 
hybrid perturbations by considering the long wavelength modes, which can be solved analytically. § [v| is devoted to 
the numerical investigation for all cosmological scales. Finally, we present our conclusions in §|vj| There is also an 
appendix, which details the evolution equations and gives the full solution for the long wavelength modes. 

II. ADIABATIC AND ISOCURVATURE PERTURBATIONS 

In this section, we will define precisely the notion of adiabatic and isocurvature perturbations, and introduce the 
notation that will be used throughout this paper. Only perturbations of the scalar type will be considered here. 

The perturbations of the spacetime geometry will be described by two scalar potentials, <& and which appear in 
the linear perturbation of the (flat) Friedmann-Lemaitre-Robertson- Walker (FLRW) metric, 

ds 2 = -a 2 (?7)(l + 2$)d?y 2 + a 2 (?7)(l - 2^)S i:j dx l dx j , (1) 

choice which corresponds to the longitudinal gauge (77 is the conformal time and is the Kronecker symbol). 

For matter, we will consider four different species : two relativistic species, photons and massless neutrinos; two non- 
relativistic species, baryons and cold dark matter (CDM). Their respective energy density contrasts will be denoted 
5 7 , 5„, 5b and S c (8a = Spa/pa)- 

Before defining adiabatic and isocurvature perturbations, let us introduce the idea of entropy perturbation, which 
can be defined for any pair of components, A and B, by the expression 

_ 5n A 5n B , . 

£>a,b = , (2) 

n A ub 

where ha represents the particle number density for the species A. When the equation of state for a given species is 
such that w = p/ p = Const, then one can reexpress the entropy perturbation in terms of the density contrast, in the 
form 

Sa,b ^ - T ^-. (3) 

1 + WA 1 + Wb 

By definition, an adiabatic (or iscntropic) perturbation corresponds to the case where all the entropy perturbations 
are zero. For our four components, this can be expressed as : 

3 3 

&c = h = -8 U = -<5 7 = 8 {a) , (4) 

where <5( a ) will denote the common value (up to the adiabatic index of the equation of state) of the density contrasts. 
An adiabatic perturbation is thus characterized by a unique amplitude, which can be (5( a ) but which is usually given, 
for convenience, in terms of the gravitational potential 'J', which can be directly relation to <5( a ) through Poisson's 
equation (see appendix). 

An isocurvature perturbation, as indicated by its name, corresponds to a perturbation for which the gravitational 
potential perturbation is zero (or approximately zero). To get a non-trivial isocurvature perturbation, one must 
therefore have several components and at least one non- vanishing entropy perturbation Sa,b- For N species, there 
will be one adiabatic mode and N — 1 independent isocurvature modes (there is also the possibility to have velocity- 
type isocurvature modes, see |Q). It is to be noticed that the adiabatic or isocurvature character of perturbations 
is not time-invariant. In the cosmological context, when one talks about adiabatic or isocurvature perturbations, one 



2 



implicitly assumes that this property corresponds to the initial state of the perturbations, which means deep in the 
radiation era when the wavelength of the perturbation is much larger than the Hubble radius. 

In the following, we shall consider primordial perturbations which are hybrid, i.e. which are a sum of adiabatic and 
isocurvature modes. In order to simplify the exploration of the parameter space, we will restrict our attention to the 
case of perturbations for which all entropy perturbations within three of the species are zero whereas the last species, 
X say, departs from adiabaticity. This means that three of the four species will satisfy the above relation (|j) while 
the contrast density for the remaining species X will be written in the form 

S x 

= S X +S(a), (5) 



1 + WX 



expression which defines the isocurvature perturbation Sx associated with the species X. Varying X, one can 
construct four hybrid perturbations of this type. 



III. CORRELATIONS 



In cosmology, perturbations are treated as homogeneous and isotropic random fields. It is convenient to deal with 
them in Fourier space rather than ordinary space, and all quantities defined previously can be transformed into their 
Fourier components, according to the relation (since we work only in flat space) : 

Wo^ !kx/(x) - (6) 

Primordial perturbations are usually assumed to be Gaussian, in which case their statistical properties can be 
summarized simply in terms of their power spectrum, defined for a quantity / by 

(/ k /^)=2vr 2 fc- 3 P / (fc) ( 5(k-k'). (7) 

When primordial perturbations are described by several quantities, such as would be the case if one has a mixture 
of adiabatic and isocurvature perturbations, one can also define, for any pair of random fields / and g, a covariance 
spectrum Cf >g (k) by the following expression : 

= 2^fc- 3 C/, 9 (fc)<5(k - k'). (8) 

Correlation between / and g can also be expressed in terms of a correlation spectrum Cf g {k) obtained by normalizing 
Cf,,(k) : 

In the present work, we are especially interested in the possible correlation between adiabatic and isocurvature 
primordial perturbations, i.e. C$,s(fc). Until very recently, only independent mixtures, i.e. with vanishing correlation, 
were considered in the literature. This statistical independence means that the quantities <!> and S can be expressed 
as 

<f = Vl /2 e l7 S = vl /2 e 2 , (10) 

where e\ and ei are independent normalized centered Gaussian random fields (i.e. such that (e 2 ;(k)) = 0, 
(ei(k)e*(k')) = <5y £(k — k'), for i, j — 1, 2), and where the subscripts k are implicit, as it will be the case in the rest 
of this paper. With the assumption (|Io|), one obtains immediately vanishing covariance and correlation spectra. 

However, as it was shown in a specific model of double inflation [jul, one can also envisage models of the early 
universe where correlated primordial perturbations are generated. To be more specific, this would be the case if 
one imagines several independent stochastic processes taking place in the early Universe, which contribute both to 
adiabatic and isocurvature perturbations, i.e. such that 



$ = ]T$,; e ,, S = ^2Sie h (11) 
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where the are independent normalized centered Gaussian random fields. In the specific example of |10| , there were 
two independent random fields, generated by the quantum fluctuations of two scalar fields. 

In the present article, our goal will be to study systematically the consequences of a totally correlated mixture of 
adiabatic and isocurvature perturbations, i.e. primordial perturbations which can be written in terms of one single 
random field. Of course, the consequences of more general initial conditions can then be obtained by simply adding 
the spectra (to get the total density power spectrum or the total temperature fluctuation multipole spectrum) of 
several totally correlated initial conditions. 



IV. LONG WAVELENGTH ANALYSIS 



As shown in the appendix, it is possible to solve analytically for the long wavelength perturbations. Totally 
correlated perturbations can be defined by two primordial quantities : the gravitational potential perturbation deep 
in the radiation era, denoted $, and the entropy perturbation denoted Sx (the index X depends on the species which 
departs from adiabaticity as explained in Section ||). It is then possible, using the expressions of the appendix, to 
compute observational quantities at the time of last scattering as functions of the two primordial quantities <J> and 
Sx , for scales larger than the Hubble radius at the time of last scattering. 

In this section, we will use, instead of the energy density contrasts 5a defined previously in the longitudinal gauge 
(IT]), the slightly redefined energy density contrasts (defined in the flat-slicing gauge) : 



4*, 



5 V - 4*, A b 



3*, 



3*. 



(12) 



The reason to use these quantities is essentially that the conservation equations look much simpler (see appendix). 
Moreover, it is to be noticed that, with these new density contrasts, the definitions of adiabatic and isocurvature 
fluctuations keep exactly the same form. In other words, a purely adiabatic perturbation, as defined by (^), will also 
be characterized by : 



-A, 



Ho)- 



(13) 



For a mixed perturbation, with the species X deviating from adiabaticity, the same relation will hold for the three 
species other than X, and the density of the latter will be given by 



1 + Wx 



S 



X 



Ha)- 



(14) 



Our purpose will now be to express the observable quantities for long wavelength modes, namely the gravitational 
potential in the matter era and the temperature anisotropics, in terms of the primordial quantities $ and Sx- This 
will be possible by using the two following relations, which are demonstrated in the appendix. The first relation gives 
the metric perturbation in terms of the primordial density contrast, during the radiation era, 



ad 



4> 



2 1 
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the second one being the analogous equation during the matter era, 

1 

5 



^matter = ~ r (^b^Afc 



^ D A C ) 



(15) 



(16) 



Note that, in the above equations, f2^ D and Sl!^ D are taken in the radiation era, whereas O^ 10 and f2^ D correspond 
to their values in the matter era. 

Our hybrid perturbations can be specified either by the pair ($, Sx) or the pair (A/ a \, Sx)i the relation between 
the two following immediately from the relation ([l5|). For a purely adiabatic perturbation, it is easy to see, using 
(H), that 



A(a) = -- I 3 



5 " 



(17) 



For a mixed perturbation, there will be in general an additional term proportional to Sx ■ I n the case of the baryons 
and CDM, the expression for A( a ) is the same as the adiabatic case, simply because and A c do not appear in (|lj 
For the relativistic species, one obtains 
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A(a) = -| (s + * - (l - ^°^7 = <** + (18) 

for a photon-type hybrid perturbation and 

A (a) = -| (3 + ^0* D ) * - J (9 - 4f£ D ) ft^S, EE a$ + (19) 

for a neutrino-type hybrid perturbation. 

Substituting in ( flq) the expressions of Af, and A c in terms of A( a ) and Sx, thus in terms of $ and Sx, it is then 
possible to find the gravitational potential perturbation during the matter era. For a purely adiabatic perturbation, 
one finds 

3 f, . 4„ RD 



$adiab = ^ (3 + -<V> *. (20) 

One recognizes the standard transfer coefficient of 9/10 if one ignores the anisotropic pressure of neutrinos (see 
e.g. |L9]]). Here, we have its generalization, which is numerically very close to 1, when the anisotropic pressure is taken 
into account. For a hybrid perturbation, the gravitational potential perturbation during the matter era, will be of 
the form 

$ = *adiab + 3>isoc, (21) 



where $ a diab corresponds to the term proportional to which is, in all cases, given by the same expression (|2C|), and 
<i>i SOC is the term proportional to Sx , whose explicit expression depends on the particular species considered. For 
baryons and CDM, its form is simply 

$isoc = -^ D 5 x , X = b,c. (22) 
5 

For a photon-type mixed perturbation, one finds 

$isoc = ^l-^ D )^ D ^ (23 ) 

and finally, for a neutrino-type hybrid perturbation, one gets 

$isoc = ^(9-4ft« D )^ D S,. (24) 

The decomposition ( pl| ) expresses the fact that a primordial isocurvature perturbation will also contribute, in the 
matter era, to the potential perturbation, whereas it is of course not the case in the radiation era. This illustrates, 
once more, that the separation between adiabatic and isocurvature modes is not conserved during time evolution. 

Let us now evalute the contribution of the primordial perturbations to the CMB temperature anisotropics, here 
only for large angular scales since we are restricted to long- wavelength perturbations. Neglecting a local monopole 
and dipole contribution, the temperature anisotropies, due to scalar perturbations, are approximatively given by : 
(see ||, 0) 



AT 1 f Vo / • • \ 

~T = tA 7LSS + (* + *)lss - e'&fliss) + / (<!> + *) dA, (25) 



where e % is a spatial unit vector corresponding to the direction of observation, the subscript LSS indicates that the 
quantities are evaluated at the last scattering surface, a dot denotes a derivation with respect to the conformal time 
77, r)o is today's conformal time, and the integral in the last term runs on the photon line-of-sight. The contribution 
due to the first two terms is usually called the Sachs- Wolfe term (SW), while the third term is called the Doppler 
term and the last one the integrated Sachs- Wolfe term (ISW). In general, but not always (see the pathological cases 
below), the SW term is dominant for large angular scales. In terms of our variables, the SW term can be written 



1 . _ , 1 



/AT\ 

V~Jsw = 4^ + $ + ^ ~ 4^ 4 ^ K l2,il 
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where the quantities are evaluated at the last scattering surface, assuming that last scattering occurs well in the 
matter era (in this case \I/ ~ 4>). Using the expressions obtained above, it is now possible to express the SW term 
as a function of the primordial perturbations and Sx- As for the gravitational potential perturbation, one can 
decompose this term into 



/ AT\ _ ( AT \ ( AT 

V T ) sw \T ) a diab V ™ 



(27) 



where the adiabatic component is the term proportional to and the isocurvature component is proportional to Sx ■ 
For all types of hybrid perturbations, the adiabatic term is the same : 

=^(3+W)4. (28) 



V t ; adiab 10 v 5 

Note that one has AT/T adiab = ^adiab/S, which reminiscent of the standard (adiabatic) Sachs- Wolfe term. As for the 
isocurvature term, it will depend on the particular type of perturbation. For hybrid perturbations which are baryon 
or CDM isocurvature, one finds 

(^T) =- 2 ~nTSx, X = h,c. (29) 

Note that, for baryons and CDM, one has the relation (AT/T)j soc = 2$ isoc , and by comparison with the similar relation 
for the adiabtic terms, one recognizes the standard statement in the literature that pure isocurvature perturbations 
(of the baryon- or CDM-type) produce large scale temperature fluctuations six times bigger than pure adiabatic 
perturbations. For photon isocurvature hybrid perturbation, one will get 



6-"-n* D + U^ D Y)s 7 , (30) 



V T J isoc 15 V 5 " 5 
whereas, for neutrino isocurvature hybrid perturbation, the expression is 

/ AT \ 1 



, (9-40^)0^. (31) 

V T /• 75 



V. NUMERICAL ANALYSIS 



The present section will be devoted to the predictions of temperature anisotropies, as well as large scale structure 
power spectrum, for primordial correlated hybrid adiabatic and isocurvature perturbations. We will keep fixed a 
certain number of parameters : Oa = 0, Oq = 1, h W Q — 0.5, fib — 0.05, three species of massless non-degenerate neu- 
trinos (leading to ^ D = 1/[1+ (21/8)(4/ll) 4 / 3 ] ~ 0.595 and = 1 - ^ D ~ 0.405), and standard recombination. 
The primordial perturbations will be assumed to be scale-invariant. 



A. Temperature anisotropies 

As far as temperature anisotropies are concerned, isocurvature perturbations can be distinguished from pure isocur- 
vature perturbations by a much larger plateau, as explained in the previous section, with the consequence that the 
all the acoustic peaks will appear smaller than this plateau. To show these two extreme behaviours we have plotted 
on Fig. the case of pure adiabatic initial conditions, and on Fig. | the case of pure CDM-type isocurvature initial 
conditions. We have also plotted, in each case, the SW, Dopplcr and ISW contributions. We have used, as usual, the 
angular power spectrum for the temperature anisotropies, defined by 

CF = (|«L| 2 ), (32) 

where the aj m are the multipole coefficients that appear in the decomposition into spherical harmonics of the tem- 
perature fluctuations, i.e. 
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(33) 



In the case of hybrid perturbations, we will be somehow between these two extreme situations, 
let us parametrize the hybrid perturbations by A, which is defined in the relation 

S = A$, 



For convenience, 



(34) 



which will quantify how far we are from a purely adiabatic model. The case A = corresponds to pure adiabatic 
initial conditions, whereas the limit where A goes to infinity corresponds to pure isocurvature perturbations. A can 
be positive or negative. To be more specific, one can call the hybrid perturbations we are studying correlated when 
A > and anticorrelated when A < 0. In Figs. |[]^, we have plotted the total temperature anisotropy as a function of 
the multipole index i, for various values of the parameter A (for CDM-type hybrid perturbations), and keeping the 
same normalization at large angular scales. 

To emphasize the difference between correlated hybrid perturbations and independent hybrid perturbations, which 
have been considered in the literature, we have plotted , in Fig. ^J, the total temperature anisotropy for independent 
hybrid initial conditions. The curves are parametrized by the number 1Z, which is defined by 



(35) 



In some sense 1Z is the analogous, in the independent case, of A since the square of both quantities corresponds to the 
ratio of the power spectra. But, of course 1Z can be only positive. The way these curves are obtained is also different. 
Whereas for the correlated mixtures, one implements hybrid initial conditions from the beginning and one runs the 
Boltzmann code once (per model), in the case of independent mixtures, one runs the code first with purely adiabatic 
initial conditions, then a second time with purely isocurvature initial conditions, and the final Ci are obtained by a 
weighted sum of the Ci obtained from each run. As a consequence, the first acoustic peak, as well as the subsequent 
ones, will always appear lower, relatively to the plateau in the hybrid case than in the purely adiabatic case. 

The behaviour of the Cg for the correlated hybrid models is quite different when one increases the isocurvature 
proportion. For anticorrelated perturbations, i.e. A < 0, the behaviour is what is expected naively : the amplitude 
of the peaks decreases, relatively to the plateau, with higher proportion of isocurvature perturbations, as illustrated 
in Fig. H with the curves lower than the adiabatic case. But the evolution is more complicated when one considers 
correlated models, i.e. with A > 0. Starting from the adiabatic case (A = 0) and increasing slowly A, one begins with 
a phase where the amplitude of the peaks will increase with respect to the plateau, as illustrated by the curves above 
the adiabatic one in Fig. [| If one goes on increasing A, one reaches a critical value, beyond which the peaks will now 
diminish with increasing A, as illustrated by the curves of Fig. |. 

One can understand this surprising behaviour if one goes back to the results of the previous section. In the case of 
CDM hybrid perturbations, one can evaluate the SW plateau, using eqns. (|2g|) and (29) : 



± / sw 



10 



-n 



RD 



-n MD x 



(36) 



and therefore, there is indeed a critical value for A for which the SW plateau is suppressed, which explains the relative 
height of the peaks. In fact, things are slightly more complicated, because when the SW term is suppressed, because 
this special choice of initial conditions, the other terms which contribute to the anisotropics cannot be neglected any 
longer. Fig. |t] shows in particular that the plateau can be due essentially to the ISW effect. Note that the value 
A = 1.36 for which this effect was numerically obtained is slightly different from the value one would deduce from (|3C|). 
This is because ( pff ) was obtained by supposing that the last scattering surface is completely in the matter-dominated 
era, which is not the case since the radiation-to-matter transition occurs no very long before recombination. One can 
also adjust the initial conditions so that the large scale anisotropies will be dominated by the Doppler term, in which 
case there is no longer a plateau on large scales, but an increasing slope as can be seen on Fig. [| 

Although at this stage we have discussed and illustrated only the CDM correlated hybrid case, a similar behaviour 
appears for the three other types of correlated hybrid perturbations, but with noticable differences. We have sys- 
tematically explored the parameter space for the four types of correlated hybrid initial conditions and measured the 
predicted height of the first acoustic peak with respect to the plateau. The results are given in Fig. |9[ Here, we have 
adopted a different parametrization of the hybrid correlated perturbations so that one can represent easily all cases. 
We have defined an angular variable 0x so that our initial conditions for the density contrasts are of the form 



^x 



l + wx 



cos Vx j 



1 + WA 



= sin 9 X , X. 



(37) 
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which implies 



S x = (cote x -l)A {a) . (38) 

Of course, this parameter Ox can be related to the parameter A. Let us write 

A (a) =a$ + p x S x , (39) 

where the subscript X for the coefficient fix refers to the type of hybrid mode we are considering, a is the same for 
all four types and is given by the coefficient in cq . (|l^ ). (3 X is zero for the CDM and baryon modes and is given by 
the second term on the right hand side in eqns. (|18|) and ( [l9| ) for the photon and neutrino cases respectively. The 
relation between Xx and 9x is then 

Xx l-(cot0 x -l)/V (40) 

In the region near — 7r/4 (corresponding to the pure adiabatic case and where all curves cross), it is easy to see that 
the relation between and A is the same for the four types of perturbations and is given numerically by 

A ~ 2a(0 - it/4) ~ 9.97(0 — 7r/4). (41) 

Let us see what happens when one deviates from the pure adiabatic case. In the baryon case, the first peak will 
increase for correlated mixtures and decrease for anticorrelated mixtures, as in the CDM case, although in a much 
more moderate way. On the contrary, in the photon and neutrino cases, the first peak will start to increase for 
anticorrelated perturbations. Moreover, the evolution in the neutrino case is slower than in the photon case. 

All these results can be understood rather easily with the analytical results of the previous section. Indeed, in all 
cases, the SW anisotropy term can be written in the form 



/ AT 



= (a + b x \x)<5>, (42) 

sw 



where the special case of CDM is given just above. The coefficient a is the same for all four cases, as was shown 
in the previous section. Therefore, the evolution of the SW plateau will be determined by the coefficient bx which 
is different in each case. For baryon and CDM perturbations, the coefficient b is negative, which explains why the 
increase of the first peak (or equivalently the decrease of the plateau) corresponds to correlated (A > 0) perturbations, 
when one deviates from the pure adiabatic case. Moreover, \b c \ > because f2 c > Of, and therefore the response 
to the increase of the isocurvature proportion is stronger in the CDM case. In the photon and neutrino cases, the 
coefficient f3 is positive and therefore only anticorrelated perturbations can lead to an increase of the peak. Moreover, 
b 7 > b v , and similarly to the heavy species, the response to the increase of A is stronger in the photon case than in 
the neutrino case. A consequence of these results is that, potentially, the correlated hybrid perturbations in the CDM 
and photon cases will be more strongly constrained by the CMB measurements than the baryon and neutrino cases. 

It is also important, in the spirit of putting constraints on this type of modes, to see the position of the two Doppler 
peaks on the £-axis. We have plotted in Fig. the positions of the first and second acoustic peaks. Near the 
pure adiabatic case, the behaviour is once more strongly pronounced in the CDM and photon cases. In this region, 
in the CDM case, correlated hybrid perturbations tend to displace the peaks to smaller I whereas anticorrelated 
perturbations push the peaks to higher I. The same behaviour, but very attenuated, seems to apply to the baryon 
and neutrino cases. In contrast, in the photon case, higher I correspond to correlated perturbations and lower I to 
anticorrelated perturbations. 

Another important feature of the acoustic peaks is the impact of the isocurvature part on the amplitude of the 
second peak. Fig. [ll] gives the evolution of the relative amplitude of the second peak with respect to the first peak, 
when one varies the coefficient A, for the four types of perturbations. As before, in the vicinity of the pure adiabatic 
point, the effects are more important in the CDM and photons cases than in the two other cases. 



B. Power spectrum 

In the previous subsection, we have analysed the CMB fluctuations. Another important part of observational data 
comes from the large-scale structures. However, it is more difficult to infer as precise information from the LSS data 
as from the expected CMB measurements, because of additional complications such as the bias effect. Moreover, the 
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signatures due to pure adiabatic perturbations and pure isocurvature perturbations respectively is easier to distinguish 
in the CMB anisotropics than in the matter power spectrum. But, nevertheless, the LSS power spectrum is a useful 
tool, at least to check the overall amplitude of the perturbations 



14 



To illustrate the power spectra corresponding to correlated hybrid perturbations, we have plotted on Figs. ^ 
the power spectra (taking the gravitational potential perturbation as reference variable) in the case of CDM-type 
correlated hybrid perturbations, for various values of A. 

Some interesting behaviour occurs for large values of A. When A > a, it is easy to see using eq. ( pio| ) that the CDM 
and the baryon density contrasts have initially opposite signs. The CDM density constrast evolves [see eqns. ( A.1 JA2|) 1 
according to the equation 



which can be solved to give 



. init 



a&dr) 



dry', 



(43) 



(44) 



where AJ. nlt is the initial value for the CDM density contrast, so that it will evolve towards values with sign opposite 
to that of that is of same sign as Ar ). Moreover shortest wavelengths will evolve rapidly enough and change sign 
whereas longest wavelengths will not, so that the CDM power spectrum should exhibit a sign change. For a given 
mode, this sign change occurs all the more rapidly as the initial CDM density contrast is small as compared to <SE>, 
that is when A is small. Plotting the matter power spectrum, it is therefore natural to expect that the wavenumbcr 
at which it is zero is all the more small as A is small. This is what we can check on Fig. [TIL 



C. CMB polarization 

Although CMB polarization has not yet been measured and is expected to be difficult to measure, it can provide 
a lot of additional information concerning the cosmological perturbations and the cosmological parameters [£2j . This 
is why we will consider briefly the consequence of correlated hybrid perturbations on the CMB polarization. The 
polarization can be decomposed | p2| into E-mode polarization and E-mode polarization, and we shall define, in 
addition to the temperature angular power spectrum, the E-mode angular power spectra, 

CF E ^(\af m f), (45) 

and the correlation spectrum 

Cf T ee (aj*af m ), (46) 

where the aj m correspond here to the same coefficients as those defined in (|3^), which gives the correlation between the 
temperature and the -E-mode polarization. Scalar-type perturbations do not contribute to the B-mode polarization. 
As an illustration, we have plotted the E-mode polarisation anisotropy spectrum for various values of A. Two features 
in these curves are obvious. First, the amplitude of the polarisation varies, as well as the position of the first peak 
in the spectrum. This can be understood by looking at the Boltzmann equation for the photon fluid. Before the last 
scattering surface, Thomson diffusion in important, and the photons anisotropic stress approximatively obeys to the 
equation : 



4 k 
15 £ 



Vby, (47) 



where k stands for the differential Thomson opacity, and VJ, 7 is the photon-baryon plasma velocity. Thus, the photons 
anisotropic stress (which is proportional to the E-type polarization, see eqns. (62,63,77) of p3]] ) is proportional 
the photon dipole, and the Cf E are proportional to the Doppler contribution of the CMB temperature anisotropy 
spectrum. Therefore, t he p osition and height of the peaks vary for both spectra in the same way with A. These results 
are represented in Fig. |l5J and the cross-correlation spectrum between temperature and polarisation can be found on 
Fig. M. 







VI. CONCLUSION 



In the present work, our goal has been to analyze the effects of correlated hybrid adiabatic and isocurvature 
perturbations on observational quantities, in particular the CMB anisotropics. We have isolated four "elementary" 
modes, corresponding to a deviation from adiabaticity of one of the four standard species (photons, neutrinos, baryons, 
CDM) : each type of these elementary modes is characterized by only two parameters (S and A( a )). One could, of 
course, generalize our analysis by relaxing adiabatic ratios among the three remaining species, but at the price of 
requiring four parameters. 

We have shown, in the case of these elementary modes, that correlation leads to very specific effects on the CMB 
anisotropies and on large scale structure which do not appear in the case of independent mixtures. 

In this paper, our purpose was not to confront directly observations with this type of models, but rather to focus 
on some qualitative interesting consequences of correlated hybrid perturbations. This is why we have considered only 
a subclass of models, which are extreme in the sense that they are totally correlated and simple because they are 
described by only two parameters. If one wishes to compare hy brid models with observations, one should consider the 
sum of spectra of the type we have obtained (as explained in § |lll| ) , which means that the adiabatic and isocurvature 
perturbations would then be only partially correlated (this is the case in the specific model of p0|). 

The present data are still too imprecise to be able to distinguish this kind of correlated hybrid perturbations, but 
it may be interesting to know how much precise data would constrain these modes. In practice, it might turn out to 
be a difficult task to disentangle the presence of such modes in the data, unless one assumes a specific early universe 
model. 
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APPENDIX A: EVOLUTION OF PERTURBATIONS 



In this appendix, we derive the evolution of all quantities for long wavelengths, i.e. scales outside the Hubble radius : 
k <C aH . The whole system of equations, in a notation slightly different from the one adopted here, can be found for 
example in [ p4| . 

Let us first introduce the system of equations governing the evolution of the matter perturbations, which related the density 
contrasts of the four species to the scalar component of their velocities (denoted by Va)- We first have four equations of 
conservation, for each of the four species, that read (in Fourier space) : 



A c 
A 7 
A 6 



W c 

-W bl 

kVb-y, 



(Al) 



where a prime denotes a derivative with respect to the conformal time and where Vt, 7 is the velocity common to the baryon 
and photon fluids, which are coupled until last scattering. We then have three Euler equations, two for the independent 
fluids of CDM and neutrinos, and one for the baryon-photon fluid : 



V u = -k 

% = -nv c 

Vtn = - 



4 



$ - o\ 



4fi 7 + 3fii 



try 



k- 



4fi 7 



4n 7 + m b 



A 7 

T 



(A2) 



where Ti is the comoving Hubble parameter defined by TC = a' /a. In the last equation, the coefficients Qa are time 
dependent since the ratio of the energy density of a given species with respect to the critical energy density, will change 
with time. 

To close the above system of equations, one needs the Einstein equations, which expresse the metric perturbations in terms 
of the matter perturbations. Only two components of the Einstein equations are useful, the other ones being redundant, 
and they are the Poisson equation 



k 9 

and the anisotropic stress equation, 



* = |5> 

X 



A x - — (l + w x )Vx 



^ (*-<*>) 



(A3) 



(A4) 



where a v represents the anisotropic stress due to the neutrinos (which thus require a decription beyond the perfect fluid 
approximation). To get the evolution of <t„, one must use a higher moment of the Boltzmann equation (see |[24j| ), 



(A5) 



where other terms on the right hand side have been neglected. 

In order to solve the above equations for long wavelengths, it is convenient to consider an expansion of all quantities in 
terms of the small parameter krj, so that one will have 



X = X <0) + X {L 'krj + X yz> {kriY + 



(2), 



(A6) 



The Euler equations (A2) then enable us to express the first order velocity term as a function of the zeroth order density 
and gravitational potentials, 



V, 



(i) 



(o) 



Co) 



(i) 



Vr 



-1a (0) 

4 " 



l+Hr, 

n 7 Ay } + 4O 7 vfr (0) + (4» 7 + 3»i,)$ (0) 
4H 7 + 3(1 + Hr))Q. b ' 



(A7) 



The zeroth order components of the velocities are as usual set to zero, otherwise one would get a divergence in the right 
hand side of (A£) (unless there is a special cancellation of the type mentioned in |lq|). 
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Using (A.4) and (|A5|) at lowest order, one gets 



5 

and susbtitutirig the above expression for one finds the following relation between $* ' and ^p' - 1 



(A8) 



i + ~n v H 2 rf ) $ (0) = ( i - in,n z v z ) * w - ~o,«va 



2 2 ,T,(0) 



V 2 ^2 A (0) 



(A9) 



Substituting in Poisson's equation (A3), at lowest order, the expressions obtained above for the velocities, and using (A9), 
one finally gets the following cumbersome equation, relating <E>' ' to the four species densities, 



(3 + 7 +0„) (l + |a,y 2 ) 

V^-a f i - t n »y 2 ) + 8 V^ \ 7 

l+y V 5 / Jl 



(4fi 7 + 3n 6 )y 



4fi 7 +3fi b (l + y) 
$(0) 



(o) 



in 7 (l + y) + 3n 6 (l + 2v) (0) 
4f2 7 + 3fi b (l + 7 7 

5 4» 7 + 3^(1 + y) + 1 + y + (3 + " 7 ~ 3 "- } y 



(0) 



(A10) 



where j/ = Tir;. While this expression yields the evolution of the gravitational potential perturbation during the whole 
evolution of the universe from the deep radiation era till the last scattering, it will be sufficient for our present purpose to 
retain from this equation only its asymptotic forms in the radiation era and the matter era. In the radiation era , y — 1 
and Q, c , f2j, <C f2 7 , Q v , so that the above expression simplifies to give 



3» 



(o) 



4 

3 + -a, 

5 



- 1 r f 4 

2 1-- 



n„) n 7 a 7 0) + 1 (9 - 4n„) a,Al 0) 



In the matter era, y — 2 and fi 7 , f2„ <C f2 c , f2{,, so that one finds 



,C°) 



(o) 



(All) 



(A12) 




FIG. 1. CMB anisotropies in the pure adiabatic model (A = 0). The solid line represents the total (scalar) contribution. The 
Sachs- Wolfe, Doppler and Integrated Sachs- Wolfe contributions are respectively represented in long-dashed, short-dashed, and 
dotted lines. At large angular scales (low £), the total amplitude is essentially due to the Sachs- Wolfe contribution. 
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FIG. 2. CMB anisotropies in the pure isocurvature CDM model (A = ±oo). The solid line represent the total scalar contri- 
bution. The Sachs- Wolfe (SW), Doppler and Integrated Sachs- Wolfe contributions are respectively represented in long-dashed, 
short-dashed, and dotted lines. Note that the power at large scales (low £) is higher than at the degree-scale. 
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FIG. 3. CMB anisotropies in CDM-type correlated hybrid models for various values of the parameter A. The highest curve is 
studied in more details in Fig. ^, the dotted curve represents the (standard) adiabatic case, and the lowest represents the pure 
isocurvature case shown in Fig. [| Note that the height of the acoustic peaks with respect to the Sachs- Wolfe plateau varies 
with A, according to eq. (k2j). 
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FIG. 4. CMB anisotropies in CDM-type correlated hybrid models for various values of the parameter A. The two highest 
curves are studied in more details in Figs. |^ and [| 
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FIG. 6. CMB anisotropies in independant (uncorrelated) hybrid CDM-type models for various values of the parameter 1Z. 
The lowest, solid curve represents the pure isocurvature case of Fig. |^ (1Z = oo), and the highest, dotted curve represents the 
adiabatic case (7Z = 0). We have chosen for 7Z the same numerical values as for A in Figs. ^J-^. Note the significant difference 
between the correlated and the independant cases especially in the region where 1Z, A ~ 1. 




FIG. 7. CMB anisotropies in a CDM-type correlated hybrid model (with A = 1.36). The solid line represents the total 
(scalar) contribution. The Sachs- Wolfe (SW), Doppler and Integrated Sachs- Wolfe contributions are respectively represented in 
long-dashed, short-dashed, and dotted lines. The parameter A has been chosen so that the SW contribution cancels at lowest 
order. In this case, the low multipole power is no longer dominated by the SW contribution, but rather by the ISW contribution. 
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FIG. 8. CMB anisotropics in a CDM-type correlated hybrid model (with A = 0.95). The solid line represent the total 
contribution. The Sachs- Wolfe (SW), Doppler and Integrated Sachs- Wolfe (ISW) contributions are respectively represented in 
long-dashed, short-dashed, and dotted lines. The parameter A has been chosen so that the SW and ISW contribution almost 
cancel each other at lowest order. In this case, the low multipole power is dominated by the Doppler contribution, which is not 
flat at low multipoles. 
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FIG. 9. Ratio of the height of the first acoustic peak to the height of the Sachs- Wolfe (SW) plateau for the four types of 
hybrid models. The height of the SW plateau is obtained by averaging the power between I — 2 and £ — 25, which correspond 
to the angular scales explored by COBE. We have defined the "first acoustic peak" as the first maximum of the multipole 
anisotropy spectrum for £ > 100. Each ratio peaks at a high value (~ 30) which roughly corresponds to the moment where the 
SW plateau disappears [see eq. (^2|); in practice, this occurs when the SW and ISW contributions cancel at low multipoles]. 
The four curves intersect at 6 — 7r/4, as expected, since this value corresponds to the pure adiabatic case [see e q. (|38[ )1. As 8 
varies, the position on the ^-axis of the first acoustic peak slightly shifts to the left or to the right (see also Fig. 10[). In some 
cases, the peak goes below £ — 100. In this case the "new" first acoustic peak position is at i ~ 200 — 300, and its height is 
different, hence the discontinuities in the curves. 
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FIG. 10. Position of the two first acoustic peaks in the four types of hybrid models. As explained in Fig. gj, the first acoustic 
peak is defined as the the first maximum of the multipole anisotropy spectrum for I > 100. In some cases {9 ~ 3 for the 
isocurvature baryon model, 6 ~ 0.5 for the isocurvature photon model and 9 ~ 1.2 for the isocurvature CDM model), the first 
peak disappears in the low multipole region, and the "new" first peak becomes the "former" second one. In some other cases, 
the peaks is smeared in the power spectrum and disappears, as it is the case for the first acoustic peak in the isocurvature 
neutrino model at 6 ~ 1.2, and at 6 ~ 1.7 for the second one. 
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FIG. 11. Relative heights of the first two acoustic peaks for the four types of hybrid models. The various discontinuities of 
the four curves originate from the discontinuities in the peak positions, explained and illustrated in Fig. [lo| 
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FIG. 12. Matter power spectra in the CDM-type hybrid models for various values of the parameter A. All the curves have 
been normalised to COBE. Note that the overall amplitude, as well as the position of the maximum varies with A. 
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FIG. 13. Matter power spectra in the CDM-type hybrid models for various values of the parameter A. All the curves have 
been normalised to COBE. 
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FIG. 14. Matter power spectra in the CDM-type hybrid models for high values of the parameter A. All the curves have been 
normalised to COBE. Note that one the power spectrum goes to zero for a critical scale, which depends on the value of A. 
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FIG. 15. CMB polarisation anisotropies in the CDM-type hybrid models for various values of the parameter A. We have 
represented the spectrum for the same values of A as in Fig. |^ (but we have omitted two of them for clarity). The amplitude 
of the spectrum decreases and the peaks shift to the right as one goes from the adiabatic case (A = 0, dotted line) to the pure 
isocurvature case (A = ±oo, solid line), as the spectrum closely follows that of the Doppler contribution to the temperature 
anisotropy spectrum (see Figs. |l| and \ 
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FIG. 16. CMB temperature-polarisation cross-correlation spectrum in the CDM-type hybrid models for various values of the 
parameter A. As inFig. [ll| we have represented the spectrum for the same values of A as in Fig. ^ and have omitted two of 
them for clarity. The amplitude of the spectrum decreases and the peaks shift to the right as one goes from the adiabatic case 
(A = 0, dotted line) to the pure isocurvature case (A = ±oo, solid line). 
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